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Abstract 

A deformation parameter of a bihamiltonian structure of hydrodynamic type is 
shown to parameterize different extensions of the AKNS hierarchy to include negative 
flows. This construction establishes a purely algebraic link between, on the one hand, 
two realizations of the first negative fiow of the AKNS model and, on the other, two- 
component generalizations of Camassa-Holm and Dym type equations. 

The two-component generalizations of Camassa-Holm and Dym type equations can 
be obtained from the negative order Hamiltonians constructed from the Lenard rela- 
tions recursively applied on the Casimir of the first Poisson bracket of hydrodynamic 
type. The positive order Hamiltonians, which follow from Lenard scheme applied on 
the Casimir of the second Poisson bracket of hydrodynamic type, are shown to coin- 
cide with the Hamiltonians of the AKNS model. The AKNS Hamiltonians give rise 
to charges conserved with respect to equations of motion of two-component Camassa- 
Holm and two-component Dym type equations. 

1 Introduction 

Recently, the celebrated shallow-water equation obtained by Camassa and Holm 

Ut — Utxx = —'iuUx + 2UxUxx + UUxxx ~ f^Ux, K, = COUSt (1) 

was extended in |2j by adding on the right hand side a term ppx with a new variable p, 
which satisfies the continuity equation pt + {up)x = 0. The model resulting from the above 
generalization first appeared in the study of deformations of the bihamiltonian structure of 
hydrodynamic type jHl E] and was coined 2-component Camassa-Holm equation. Soon after 
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its derivation the model was identified with the first negative flow of the AKNS hierarchy 

Another well-known integrable partial differential equation of interest to our study is the 
Dym-type equation El El El EDI : 

"^J^XXt ^''^-^x'^^XX I '^^'^^XXX txi'Xii^ • 

It can also be extended to two-component version by adding a term pp^ on the right hand 
side of (121). The resulting two variable system is shown here to be equivalent to the negative 
flow of one of the extensions of the AKNS model. It is also equivalent to the special limiting 
procedure of deformations of the bihamiltonian structure of hydrodynamic type. 

In this paper the following is accomplished: First, we explore the Schrodinger spectral 
problem of second order describing both 2-component Camassa-Holm and Dym type equa- 
tion for different values of the deformation parameter p. We show that this Schrodinger 
spectral problem can be cast into the linear 2x2 matrix spectral problem. Using sl(2) 
gauge invariance we transform the time evolution flow of the linear spectral problem into 
the AKNS first negative flow. We should point out that there exist several ways to extend 
the AKNS hierarchy to incorporate negative flows. These extensions are parametrized by a 
single parameter identified with p. We associate two different constructions of the negative 
flows of the AKNS hierarchy to the 2-component Camassa-Holm (for p = 1) and Dym type 
equation (for /i = 0). The result concerning the 2-component Camassa-Holm equation con- 
stitutes an algebraic version of the proof given in [2j. Using connections of the AKNS and 
deformed Sinh-Gordon models to the 2-component Camassa-Holm and Dym type equations, 
respectively, we are able to find explicit soliton solutions given in hodographic variables. 
The relation to the AKNS models allows us to construct a new chain of charges conserved 
with respect to equations of motion of two-component Camassa-Holm and two-component 
Dym type equations. For both hierarchies the modified AKNS Hamiltonians provide a tower 
of positive order Hamiltonians obtained via the underlying Lenard relations of the Poisson 
brackets of hydrodynamic type from the Casimir of the second bracket. 

In section |21 we briefly review the algebraic approach to the AKNS model and show 
how to extend the model in two different ways to negative time flows based on the zero- 
curvature identities. In section (HI we set up a class of two-component Schrodinger spectral 
problems parametrized by p. In the next section |31 we transform the Schrodinger spectral 
problem by the reciprocal transformation and linearize it. The resulting linear 2x2 matrix 
spectral problem is then transformed by an sl(2) gauge transformation into the AKNS Lax 
spectral problem. The time flows of 2-component Camassa-Holm and Dym type equation 
are shown to coincide with two different negative flows of the extended AKNS model. Our 
construction allows us to find, in subsection 14.11 explicit soliton solutions for various values 
of p. In section we reproduce equations of motion for 7^ and /i = cases in the 
setting of deformations of the bihamiltonian structure of hydrodynamic type. Remarkably, 
the Hamiltonians governing positive evolution flows of the AKNS hierarchy define conserved 
charges for the 2-component Camassa-Holm and Dym type equations. Also the conserved 
charges induced by the AKNS model satisfy among themselves the Lenard relations of the 
bihamiltonian structure of hydrodynamic type. Thus, the Hamiltonians of the bihamiltonian 
structure of hydrodynamic type connected to 2-component Camassa-Holm and Dym type 
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equations split into two chains, one of the positive order induced by the AKNS hierarchy 
and one of the negative order containing generators of the equations of motion defining both 
hierarchies. 



2 Extended AKNS model 

First, let us present the AKNS hierarchy in the setting of the sl(2) loop algebra endowed with 
homogeneous gradation defined by the operator Ad/ dA. A variable A plays a double role of 
a loop parameter of the loop algebra and a spectral parameter of the underlying hierarchy. 
The matrix Lax operator L for the AKNS hierarchy reads: 



d_ 
dy 



'A 


" 




"0 q 





-A 




r 



(3) 



where d/dy is the derivative with respect to "space" variable y. The matrix Lax operator can 

be compactly written as L = d / dy — E — Aq, with E = Xa^ and the matrix Aq = qa^ + ra-, 
where is the Pauli matrix and a± are given in terms of other Pauli matrices ai , (72 
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We work within an algebraic approach to the integrable models based on the linear spectral 
problem Ll'if) — 0, which simplifies considerably under a dressing transformation: 



d_ 
dy 



-E-Ao]e 



d_ 
dy 



-E. 



(4) 



where the dressing matrix = exp (X]i<o ^d'^^^) is an exponential in negative powers of the 
spectral parameter A on a formal loop space of sl(2). Similarly, for higher flows we obtain 



d 



n— 1 



a 



1=0 



at 



(5) 



where E^'^'> = X^a^ and terms Dn^ are obtained from projection of 
on the positive powers of A via expansion relation: 



n—l 



1=0 



These dressing relations give rise to the zero-curvature conditions for the positive flows of 
the AKNS hierarchy 



n— 1 



dy dt^ ^ 



© 



9y dt^ 



©-^ - 0. 



(6) 
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In particular, for n = 2 we obtain the second flow of the AKNS hierarchy: 

which reproduces the familiar vector non-linear Schrodinger equation. 

According to jTTj, the Hamiltonian densities of the AKNS model are defined as 

= - tr (^W A^-'^)) = - 5^ tr , (8) 

where A'^""-' are given by 

oo 
k=l 

where the symbol tr in expression (jSI) denotes a si (2) trace. We list the first two Hamiltonians. 
Inserting n = 1 in (jSJ we obtain: 

ni = - tr(E(°M(-i)) = ^ tT{Al) = rq . (9) 
Similarly, for n = 2 we obtain 

n2 = qry . (10) 

Next, we extend the AKNS model to include negative grade time evolution equations gov- 
erned by the zero-curvature equations ^2] 



|_ _ _ , - - 



0. (11) 



Here, we only consider the first negative flow with n = 1 and set for brevity s = In this 
case the compatibility equation (fTTj) reduces to 

{Ao)s - D[;^^ +[E + Ao, D^'^'^] = . (12) 

A general solution of the compatibility equation (fT^ is given by 

D^-^') = BS^-^^B-\ Ao = ByB~\ (13) 

in terms of the zero-grade group element, B, of SL(2), that satisfles equation: 

{ByB-'),= [BS^-'^B-\E] (14) 

or, equivalent ly, 

iB-'Bs)y= [£^-'\B-'EB] . (15) 

Here £^~^^ is an element of sl(2) algebra of —1 grade. 

Remarkably, the compatibility of the t_i flow with positive t„, n > 1 flows does not 
require that the matrix £^^^^ commutes with E = \a^, as pointed out in [13j and |T3j. It 
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turns out that all possible cases are parametrized by a parameter /x and fall into two main 
classes depending on whether /i takes non-zero or zero value. The corresponding generic 
choices of £^~^'^ are: 

^(-1) ^ |/^^3/4A for /i^O ^^g^ 
1 (T+ / A for /i = . 

Note that the value of determinant of 8^~^^ is equal to —jj^/lQX^ and 0, respectively. There 
exist other choices of 8^^^'^ for these values of the determinant but they only lead to the 
gauge equivalent copies of hierarchies derived from the choice fll(i|l . 



3 A class of two-component Schrodinger spectral problems 

Consider a linear system 

= (^^-mX + p'X'^^, (17) 

i^t = -{^ + u)iJx + ^Uxip , (18) 

for some arbitrary constant /i (see jHl E] for /i = 1 and [12] for fi = 0). Compatibility 
condition for the above system, yields three independent equations 

Pt = -M^, (19) 
rut = -2mux - m^u + pp^ , (20) 

IT^x p '^x "^xxx -I (21) 

corresponding to coefficients of A^, A and A° in the expansion of ipxxt ~ 4'txx = 0. Equation 
(1211) can be integrated to yield: 

m = p^u — Uxx + • (22) 

Here, k is an integration constant. For p ^ ^ that integration constant can be removed by 
transforming the system by Galilean transformation: 

' t t' t ^ ^ ^ d d 

' ' dx' dx ' dt dx' dt' 

In the primed system equation (j2m) becomes: 

vrrix' + rriti = —2mUx' — rrix'U + ppx' = —2 ipu — u^'x' + -f^j Ux' — Tn^'U + ppx' 

Next, performing a shift u ^ u — v and choosing velocity v such that v = n/2p^ eliminates 
the linear terms in and u^xx from the above equation. Clearly, the above argument works 
only for /X 7^ and from now on we put the integration constant n to zero as long as /i 7^ 0. 
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Note, that the positive constant /i^, that is different from one can be absorbed by ap- 
propriately redefining fields and derivatives. Defining A = A//i^, p = p/i, u = uiJ? and new 
variables x and t such that dx = fJ^dx, dt = (l//i)c^ allows us to rewrite a linear system 
(HHl-dlHl) as 

^xx = (^^-mX + p'X'^^^, (23) 

if^i = -{^ + u)^x + l:Ux^p- (24) 
2A ^ 

with m = /i^M — Uxx + = u — Uxx + ^/2. Thus, for /i^ 7^ 1 the spectral system has been 
transformed to the canonical system with /x^ = 1. 

In case of a negative /i^ (imaginary /i) we make the changes as above but with instead 
of /i^ and arrive at 

fe= (^-i-mA + p2A2^^. (25) 

Thus, only three cases of /x^ = —1, 0, 1 need to be considered separately, as concerns equations 
of motion. 

The case of /i^ = 1 corresponds to the 2-component Camassa-Holm model with m = 
u — Uxx + '^/2, introduced in while the case of /i^ = —1 corresponds to equation (pUj) 

with m = —u — Uxx + '«/2, which for p = is well-known to possess the compacton solutions 

m- 

For p = 0, we obtain from (j22p m = —Uxx + n/2. Inserting this value of m into equation 
(EOl) yields 

Uxxt = '^UxUxx UUxxx + f^Ux PPx ■ (26) 

For p = this is the Dym type equation ((21). 

After one integration (and ignoring the integration constant) we obtain from (j26p 

1 2 1 2 
=Uxt + UUxx - KU+ -U^ + -p 

^ (27) 

= {Ut + UUx)^ -KU+- [-ul + p^) . 

In terms of a new function 

^ = \ {-< + (28) 
we can cast equation ^I7\ in the following form 

{Ut + UUx)^ - KU + V = Q . (29) 

In addition, it follows from equations ^I7\ and (fTTHl that v defined by relation also 
satisfies 

Vt + {u{v + UK/2)).^ = Q, (30) 
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which becomes a continuity equation in the k = hmit. The hnear system corresponding to 
equations and (jHOJ takes a form 



(31) 



4 Transformation to the first order spectral problem. Algebraic Con- 
nection to the AKNS model 



Now, for an arbitrary /x we perform a reciprocal transformation {x, t) i-^ {y, s) defined by 
relations 

dy = p dx — pu dt, ds = dt, (32) 

and 

(9 fl fl fl fl 

(33) 



ds ^ dy 



d d d 
dx ^ dy' dt 

The commutativity of derivatives with respect to s and y variables is ensured by the con- 
tinuity equation (fT^ . Applying the reciprocal transformation and then redefining if) by 
if = ^/p4' as in [2j leads from the spectral problem (fT7j) - (fT8|) to : 



where 



fyy 
P 

Q 



(A^-PA-g) if, 

P , Py 
-2A^^ + 4A^' 



(34) 
(35) 



m 
7^ 



2 2 

_P Pyy _|_ Py 

4p2 2p 4p2 



(36) 



Our main point in this section is that we can rewrite the second-order spectral problem 
first-order linear problem: 



(37) 
(38) 







= A 






y 




Jl. 






= D 






s 







involving sl(2) matrices: 
A = 



9 

X-P 



A 



^P^i, 



Xai + gas - Pa 
Do-- 



^as + ^iPp-2g,)a.. 



(39) 
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Note that determinant of Dq is equal to detDo = — and, therefore, the matrix Dq 
becomes singular for /i = 0. Eliminating rj from the linear spectral problem (j?fj) -()38 p repro- 
duces equations (j34j) - (j35|) for ip providing that function gijj, s) appearing in (j39|l satisfies the 
Riccati equation 

Q = -g"^ - 9y 

for Q given in equation ()36|). Remarkably, the solution to the above Riccati equation takes 
a local form of 

I . (40) 

The zero-curvature equation 

As-Dy + [A, D] = 0, (41) 
can easily be derived from the linear spectral problem (|n7 j) -(|nH |) . It is equivalent to equations: 

Ps = Py, Qs + \PyP + Ppy = . (42) 

These equations were found in [2] directly from compatibility of equations It 
follows from the first of the above equations that there exists a function f{y, s) such that 
P = fy and p = fs. 

By plugging p = fs and P = fy into the second relation in ()42|1 one obtains as in j2j: 

2 fss r f \ }_f f _ fssyy I fssyfsy fssfsyy _ fssf sy _ „ / .„x 

2/^ 2 ^ 2/ 2/^ 2/2 2f^ — u. {^o) 

This appears to be the only condition, which the function / has to satisfy in order to be a 
solution of the model. 

For n = equation (j43j) simplifies to 

f^fy-fssy + ^) =0. (44) 

Is / y 

Integrating the above equation once and setting the integration constant to k/2 (see expla- 
nation below) yields 



fssy fssfsy .2^ 
J s J s J s 



+ ^ = fy, (45) 

Js Js 

or 



i^^fs)sy + \'^/fs = fsfy 

Indeed, multiplying both sides of equation by f"^ and taking a derivative with respect 
to y yields ()44|) . It remains to be shown that the choice of k/2 as the integration constant 
in equation ()45|) was consistent with equations of motion. To do this we start by recalling 
that P = m/ with 

m = -U:rx + = -P{PUy)y + o'^ 
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in the fi = case. The continuity equation ()19|) reads in the hodographic variables ps = 
—p^Uy. Accordingly, substituting Uy = (l/p)s into P we get 

D_ ips/p)y , _ r 

which is precisely equation (jl5|) . 

Let us turn our attention back to the zero-curvature equation (jH]). This equation is 
invariant under the sl(2) gauge transformation: 

A — >UAU-^ + UyU-\ D — >UDU-'^ + UsU-K 

This invariance will be used in what follows to cast the linear spectral problem ()37|) - ()38p in 
the standard form of the first positive and first negative flow equations of the sl(2) AKNS 
hierarchy. 

As a first step we gauge away the term — |po"i of order A" in the expression for D in 
equation (jH^ by choosing 

fl\f f 
U = exp I -f{y, s) 0"! = cosh — + cti sinh ■ 



X ' 7 2 2 

so that UgU^^ — \po'i = 0, due to fs = p. Consequently, 

A UAU-^ + UyU-'^ = U (AcTi + gas - Pa^) U''^ + ^Pai 

= \ai + as (^g cosh f - sinh - i (Ts (^9 sinh f - cosh / 

D —U DqU^^ = — [(— o"i cosh / — i (72 sinh /) p 
A 4A 

+ (Pp - 2gs) (as + ai cosh / + i a2 sinh /)] 

Note that the gauge transformed of the matrix D is now proportional to 1/A. 

Next, we define the constant matrix ^ = + crs), that by a similarity transformation 

maps ai to as, QaiQ~^ = —as. Note also that Qa2fl~^ = —<J2- The combined gauge 
transformations first by U and then by Q produce the final result 

A^ E + Ao = n [UAU-^ + UyU-^] n-'^ 

= Aas + ai f 5( cosh f - sinh / j + i a2 (g sinh f - ^P cosh / 



A 4A 



(46) 

iPp-2gs)as 



+ai \^{Pp - 2gs) sinh f - p cosh fj + i a2 \^{Pp - 2gs) cosh f - p sinh / 

In the above, we re- introduced E = Aas and Aq = ra^ + ga^. Comparing with the right 
hand side of equation ()46p we find that 

? = eM^?-^)' r = e-ng+^). (47) 



2 V 2 



Furthermore, defining matrix entries of as 



A 



a (3 
7 —a 



(4J 



we find from 



tliat a, (3 and 7 are given by 



(49) 



7 



-e I a + — 



Tliey satisfy tlie determinant formula + /37 = — /i^/16. Tlie matrix entries of Aq and 
enter tfie following simple identities 



2a = (3e ^ — 76-^, — ^ = /5e + 76 



/ 



P = re^ — qe 



-f 



g = - (re^ + ge ^) . 



(50) 
(51) 



It follows that the linear spectral problem (|n7 j) -(|nH j) has been transformed by the above gauge 
transformation to: 



A 



A 
-A 

a (3 
7 —a 



^ + 



q 
r 



(52) 
(53) 



for some two-component object 

We recognize in (|52j) the spectral problem L (\1') = with the AKNS Lax operator given 
by equation (jS)). It also follows easily that equation (|12j) is the compatibility equation of the 
spectral equations (j^ - lj^Hj) . The compatibility equation flT^ yields 



qs = -2/5, = 27 . 



(54) 



when projected on zero grade, and 



Py = —2aq 
'jy = 2ar , 



(55) 



when projected on —1 grade. Equations can also be directly derived from equations of 
motion 
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4.1 Examples and Solutions 



Let us recall that a general solution of the compatibility equation (fT2|) is given by expressions 
from equation It is convenient to parametrize the SL(2) group element B appearing in 
expressions in fllHj) by the sl(2) algebra elements through the Gauss decomposition: 

B = e^^"- e^^ 6^"+ . (56) 



4.1.1 The case of ^ 

As an example, we first consider jj?' = A with 8^~^'^ = (7^/2\ according to equation pfij) . As 
in ^21 in order to match the number of independent modes in the matrix Aq we impose two 
"diagonal" constraints Tr (ByB'^cr^) = and Tr {B^^Bga^) = 0, which effectively eliminate 
R in terms of and x- 

From ByB^^ = qa+ + rcr_ we obtain the following representation for q and r: 

q = -^e^ ; r = -hye-^ (57) 

where 

h = <pe^ ; h = xe^ ; A = 1 + /i /i (58) 
with a non-local field R being determined in terms h and h from the "diagonal" constraints: 

hh 

Tr (5,5- V3) = Q ^ Ry = -^, (59) 



A ' 
A~ 



hh 

Tr (5-15,^3) = Q^R, = -±. (60) 



The zero curvature equations are in this parameterization 



qs = (-%e^) =-2/^e«, (61) 



A 

r, = [-hye'^)^ = -2h^e-^ . (62) 

The two-parameter solution to the above equations can be deduced from a method combining 
dressing and vertex techniques as described e.g. in [T7j. The explicit expression is found to 
be given by: 

feexp (s/72 - ?/72) 



l + rexp(s (^-J-) -y(72-7i)) 
aexp (-s/71 + y7i) 

1 + r exp - - y (72 - 7i)) (^3) 

1 + yexp (. J-)- ^(72-70) 

1 + r exp (^s (^:^ - - y (72 - 71) 
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where 

^ _ 067172 
(71 - 72)^ 

is given in terms of four arbitrary constants a, b, 71, 72- Higher muhi-sohton solutions can be 
obtained using the same straightforward procedure. 
Comparing with equations ()54|) we find that 

he^ = e^a_, hAe~^ = a+ (64) 

where we introduced the notation 

1 

a+ = a ± — . 

By multiplying the above two relations we find that 

hJiA = (A - 1) A = a+a- = _ 1 _ 

Solutions to this quadratic equation are 

A = ±a±, hh = A-l = ±a^ . (65) 
Adding two relations in ()64j) we get 

2a = a+ + a_ = /ix H — JiA, x = e^~^ . 

X 

Solving this quadratic equation yields: 

/± = i? + ln/i-lna±. (66) 

Equation contains two solutions. The first one, namely, a+ = A, a_ = A — 1, when 
inserted into equation (jUH)) yields 

/+ = i? + In /i - In A, /_ = i? + In /i - ln(A - I) = R-\nh. 

while the second solution, = —(A — 1), = —A, leads to 

/+ = i? + ln/i-ln(-A + l), /_ = i? + ln/i-ln(-A). 

Thus, the Backlund transformation 

/+ = /_ + eln(^^^^, e = ±l 

relates the two values /+ and /_. In the reduced case of sinh-Gordon equation with h = h 
we find that A = 1 + h'^ and 

1 1 
= -In (1 + /i^) = -InA, ln(A - 1) = 21n/i. 

It follows that all the above values of /+, /_ can be summarized as: 

/. = 2ln^, e = ±l. 
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4.1.2 The case of /i = 

In the case of yU = 0, the matrix D^^^^ from equation (j^8|) takes a simple form 



a 
A 



3-/ _1 



which according to definition is reproduced by 



X 



-X 1' 
'X^ X 



for 



X 



or R= if + In a) /2. 

From equation we find for fi = that Oy = —2ag. Therefore 

q = = (InQ;)^ 

and by comparing with definition (|inj) we conclude that 

(p«), = o. 

Next, we calculate 

By = {Ry - x0,e2^) da + 0,6^^^+ + {xy + 2xRy - x'^j/e^^) a. 
Imposing condition Tr {ByB^^a^) = implies Ry — x4'y^^^ = or 

^. = -Va = -^(/. + ^)4(j-/>/2). 
What remains of expression (jTlj) is now given by 

By B-^ = (t)ye^^a+ + {Xy + XRy) ^- = {Ry/x)(^+ + {Xy + XRy) f^- 

RecaUing relations (jUHj) and (jUU]) we obtain the desired results 



that reproduce expressions (jUj). 
The compatibility equation 
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-2 
-2x' 



yields the following equations of motion 

^] = -2e'^ (72) 

X J s 

Xys + 2XsRy = ^ {xse'^)^ = 0. (73) 
Equation (|73|l implies that 

= cs{s) e-'^ , (74) 

where C3(s) is a an arbitrary function of s only. 
From equations ()68|) we find that 

Xs = fse'^ = pe'^ = C3{s)a^^e'^ 

and therefore p = 03(5) /a in agreement with relation (fTUj) . It follows that c^^s) has to be 
different from zero for consistency of the model with p 7^ 0. 
Integrating relations (f7i|) and (f7^ leads to 



X 
R 



l^sis) e-'^ ds + C2{y) (75) 



^ = -2 I ds + c^{y) , (76) 

where ci, C2 depend at most on y. Combining these two equations and setting €3(3) to be a 
constant C3 we get the deformed sinh-Gordon equation for R [18] : 

/>S PS 



Ry = -2c3 y e~^" ds J e^" ds - 2c2 / e^"^ ds + C1C3 / e"^"^ ds + C1C2 (77) 

or 

-.-2R J „ I „2R j„\ n„ „2R I „ „ -2R 



Ry^ = -2c3 e"^^ ds J ds j - 2c2e^" + 01036" 
The one soliton solution to the above equation with ci = C2 = is given by (see also [THj ) 

Ris, y) = - + 203, + In (^ ^^^.,.^.,, _ ^ , J , (78) 
where ^ 

^2 = 8C3 V 

Ki — 1 

and ko, ki are real constants. The corresponding expression for x is 

The above function x together with R from (|7H|l solve equations (f7^ - (f7!?jl . The function / 
defined by equation (f?^ provides a one-soliton solution to equation (jUj). It satisfies 

= ^ 

a 
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with 

({ki - ly/ih + 1) + koe^^'+'^^y) (koe^^'+'''y + ki + l) 

a = — — 5 -C3 . (80) 

(^koe''^^+k2y - + if 

Eliminating a from equation (f?^ using (jHE)) we get 

4C3 

/. ■ 

Therefore, comparing with we see that k = 8C3 and / given in (f7^ satisfies ()45|) and 
therefore also (HH). One-sohton solution R{s,y) given by expression (fTHj) satisfies therefore: 

K, = -- I e-'^ ds I ds. 



- (In + = 4a 



5 Bihamiltonian structure 

5.1 Bihamiltonian structure of the 2-component Camassa-Holm 
model 

As in P|, we consider the following bihamiltonian structure: 

{wi{x),Wi{x')}i = {w2{x),W2{x')}i = 0, 



{wi(x),W2(x')}i = 6'(x — x) 5"{x — x). 



(81) 



{wi{x),Wi{x')}2 = 25' {x — x), 
{wi{x), W2{x'))}2 = Wi{x)6'{x — x') + w[{x)6{x — x'), 
{w2{x),W2{x')}2 = W2{x)5'{x - x') + {w2{x)5{x - x')) , 

where now plays a role of the deformation parameter. There exists an hierarchy of 
Hamiltonians related through Lenard relations [31 

K(x), H.^}2 = jK(x), J = 1, 2, 3, ... . (82) 

The flows of the bihamiltonian hierarchy are generated by the Hamiltonians H^j^j < via: 

—^ = {w^ix),H^,h, j = 3,..., ^ = 1,2. (83) 

The lower Hamiltonians H_j for j > 1 can be obtained recursively from the Casimir H_i = 
J W2{x)dx of the first Poisson bracket applying the Lenard relations (jS^ . Following jSlHI, 
we introduce objects ipi, ip2 defined by wi = ipi — ipi^^/fJ', W2 = (p2 + V^2,x/At- Then 

{(pi{x),W2{x')}i = 6'{x - x'), {wi{x),ip2{x')}i = S'{x - x') 

and the Lenard relations yield 

H-2 = j[(f2{^i- ^i,x/lJ')]dx 



H-3 = ^ y ^2 + '^2^1 (V^i - V5i,x//u)] da; . 
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Plugging the above if_3 into equation of motion (jHHj) for j = 3 we obtain: 



{W2)t = ( ^1^2 + ^V^lV'S 



1 ^ 



X 

,2/, ,2 I „,,2 



where t = Defining u such that (pi = 2u and p such that = —p /fi + wj/A or 
= wIp'^/A — W2p'^ we can rewrite the above system of equations after a transformation 
t — > — t as 

ut-u^xt/p^ = ppx/p^ -^uUx + 2u^Uxx/p^ + uu^^^/p^ (85) 
pt = -{up)^, (86) 

which agrees with the 2-component Camassa-Holm equation. Multiplying equation (j85j) by 
p"^ and taking /i^ — *■ yields 

~Uxxt = PPx + '^UxUxx + UUxxx (87) 

corresponding to eq. with k = 0. 

In order to take the p ^ limit of the Poisson structure (|HT|) it is convenient to change 
the variables from wi, W2 to m and p defined as 

m = ^p^{wi{x) +Wi^x/fJ') = P^u-Uxx ^gg^ 
p2 = ^2 ^^2^4 _ ^ 

In terms of m and p the Poisson bracket structure (jHT|) turns into: 
{m(x), m(x')}i = 0, 

{p^{x), p^{x')}i = —p^{2m{x)5'{x — y) + m'{x)S{x — y)), 

{p2(x), m(x')}i = ^p^ {-pH\x - x') + 6"'{x - x')) , 

1 (89) 

{m{x),m{x')}2 = -p^ {p^5'{x — x) — 5"'{x — x')) , 

{p(a;),p(x'))}2 = -]^p^6'{x - x'), 
{p{x),m{x')}2 = 0, 

5.2 The fi ^ limit and the Dym type hierarchy. 

Redefining the brackets as follows 
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and taking /i — >■ limit in equation ()89|) we find for tlie first and second bracket structure in 
terms of u and p (see also [131 flQij): 

{u{x), u{x')}i = 0, 
{p^(x), p^(x')}i = 2u"{x')5'{x — x') — u"{x')5{x — x) 

{p^{x),u{x')}i = -^S'{x - x'), 

1 , , (90) 

{u{x),u{x)}2 = --d^ 6{x-x), 

{p{x),p{x')}2 = -^dj{x - x'), 

{u{x),p{x')}2 = 0, 

The first bracket in (j9(J|) has the Casimir: 



H^-l = j[p\x)-ul{x)]dx. 



This Casimir leads via Lenard relation (1821) to the Hamiltonian: 



H_2 = 2 jip'^ -ul)udx 



which in turn generates equations of motion of p = case via equations (|Hn|) and (|H^: 

^ = 7;{uxix),H.2}2 = —{uxf - uu^^ - ]-p^ 
dp I (91) 
^ = 2^Pix^,H^2}2 = - (up)^ 

This Hamiltonian structure can be extended by an additional term: 

H^2 = -2k J u^dx. 

Adding this term to H^2 will lead via relations ()91|) to correct equations of motion ()29 |) - ()30|) . 
5.3 Hamiltonians of positive order 

There exists another class of conserved charges, different from the chain of Hamiltonians 
H^j, j = 1,2, ... of negative order discussed above. These are the Hamiltonians of positive 
order originating from the Casimir 

H^^l = 2 J p{x) dx (92) 



of the second Poisson bracket (j89jl . We now employ Lenard relations (j82j) to construct higher 
order Hamiltonians. The first recurrence step: 



{■, j p{x) dx}i = {■,Ho}2 
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where "■" stands for phase space variables m{x) and p{x), leads to a new Hamiltonian: 



Ho = -r^dx. (93) 

in agreement with expression found in j4j. The integrand of Hq can be rewritten as 

mm 

— = P— = pP = Ply = fx 
P P^ 

and therefore Hq appears to be a surface term that would vanish if / would be a local field. 
On the next level we find from the Lenard relations: 

{■5 Ho}i = {-, Hi}2 , 

with 

This recurrence process can be continued to yield higher order Hamiltonians. Technical 
calculations involved in obtaining higher order Hamiltonians become increasingly tedious. 
Remarkably, we can bypass these difficulties by relying on the underlying AKNS structure 
governing higher positive flows. We recall the Hamiltonian densities Tin ((HI) of the AKNS 
model generating the positive flows of the model. Their conservation law with respect to the 
negative flow s takes a form 

(^n), = Xy , (95) 

where s and y are "reciprocal" variables describing time and space of the AKNS model and 
X is some local quantity. The above relation ensures that (J 7i„dy)s = (for X local in u 
and p) and thus the integral / Hn dy is conserved. 

In terms of the original t, x variables the conservation laws (j95|) read 

— + pu— rin = TT-X 



dt dy J dy 
or 

d ^, (9^, d d 

— Hn = -pU—Hn + ^X = -U—Hn + '^X , 

ot oy oy ox p ox 

where we used that d/dy = pd/dx. It follows that 

d d 

— {pHn) = -M^Hn - UpHnx + X^ = — {X - UpHr,) . (96) 

ot Ox 

Thus, the quantities pTin are conserved charges of the 2-component Camassa-Holm and 
2-component Dym type models. 

The first two Hamiltonian densities rq and rqy of the AKNS model, given by relations 
(jni) and (frnjl . give rise, after use of definitions (jUj) , (jHEl) and pUjl . to the following conserved 
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charges: 




which we have verified exphcitly to be conserved under equations of motion (jl9|) - (j2(J|) for 
m = fi'^u — Uxx for /i 7^ and m = —u^x + ^1'^ for /i = 0. 

We recognize in pTii the Hamiltonian Hi derived in ()94j) from the Casimir of the second 
bracket via Lenard recursion relations. Furthermore, we have shown that Hi = J pTii dx 
and H2 = J p7^2dx/2 also are interrelated via the Lenard relation: 

{■,Hi}i = {-,H2}2. 

Therefore, the conclusion is that the AKNS induced Hamiltonians pTin form the sequence of 
positive order Hamiltonians of the 2-component Camassa-Holm and 2-component Dym type 
hierarchies. The formula (jHI) given in section |21 can be used to systematically derive all the 
Hamiltonians governing positive fiows of this model. 
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